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ABSTRACT 
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Observations of Fe XVIII and Fe XIX X-ray, EUV, and FUV line emission, 
formed at the peak of Capella’s {a Aurigae) emission measure distribution and 
ubiquitous in spectra of many cool stars and galaxies, provide a unique oppor¬ 
tunity to test the robustness of Fe XVIII and Fe XIX spectral models. The 
Astrophysical Plasma Emission Code (APEC) is used to identify over 35 lines 
from these two ions alone, and to compare model predictions with spectra ob¬ 
tained with the Chandra Low Energy Transmission Grating and High Energy 
Transmission Grating Spectrometers, the Far Ultraviolet Spectroscopic Explorer, 
and the Extreme Ultraviolet Explorer. Some flux discrepancies larger than factors 
of two are found between observations of Fe XVIII and Fe XIX lines and predic¬ 
tions by APEC and other models in common usage. In particular the X-ray 
resonance lines for both ions are stronger than predicted by all models relative 
to the EUV resonance lines. The multiwavelength observations demonstrate the 
importance of including dielectronic recombination and proton impact excita¬ 
tion, and of using accurate wavelengths in spectral codes. These ions provide 
important diagnostic tools for 10^ K plasmas currently observed with Chandra, 
XMM-Newton, and FUSE. 

Subject headings: atomic data — atomic processes — stars: individual (Capella) 
— ultraviolet: stars — X-rays: stars 


1. Introduction 

The Capella system (HD 34029, a Aurigae), consisting principally of two cool giant stars 
(G8 HI -|- G1 HI), is one of the strongest coronal X-ray sources, and offers an opportunity 
to benchmark the models used in the interpretation of X-ray spectra from astrophysical 
plasmas. Although plasma codes have gone through major improvements (e.g., Mewe et al. 



1995; Brickhouse et al. 1995; Smith et al. 2001; Young et al. 2003), their accuracy and 
completeness for diagnostic analysis at high spectral resolution has yet to be fully assessed. 
Brickhouse & Drake (2000) proposed the use of the Chandra grating calibration observations 
for a series of studies, collectively known as the Emission Line Project (ELP), to test spectral 
models of collisionally ionized plasmas. This paper compares current spectral models for 
Fe XVIII and Fe XIX with Chandra observations of Capella as a step towards ensuring that 
astrophysical interpretations of spectra are based on a sound understanding of the physical 
processes involved. The Capella spectrum is well studied, showing no evidence for flares 
(Brinkman et al. 2001; Canizares et al. 2000; Brickhouse et al. 2000). The emission 
measure distribution (EMD) of the Capella system shows a strong narrow peak at 6 MK, 
near the temperature of peak emissivity for Fe XVIII and Fe XIX, producing numerous strong 
transitions. Furthermore, line fluxes from these ions show only modest variability (~20%) 
over timescales of months to years (Dupree et al. 2005), validating the combined analysis of 
multiple observations with Chandra. 

Although Fe L-shell (i.e. Fe XVII to Fe XXIV) X-ray lines offer powerful diagnos¬ 
tic potential for collisionally ionized plasmas, two long-standing atomic modeling problems 
originate from solar observations of neon-like Fe XVII (see Saba et al. 1999) that have only 
recently been addressed by laboratory programs (Brown et al. 1998; Laming et al. 2000; 
Beiersdorfer et al. 2002). Similar discrepancies between models and observations of other 
ions are now arising from Chandra spectra. For example, Xu et al. (2002) End that the 
observed Fe XVIII 3s-2p/3d-2p ratio in the elliptical galaxy NGC 4636 is higher than pre¬ 
dicted by APEC and similar to the ratio observed in Capella. The analogous Fe XVII ratio 
shows the same pattern of discrepancy, suggesting a common atomic physics origin. The 
ELP observations of Capella with three instruments offer a unique opportunity to compare 
models and observations over a broad spectral range. Comparisons among FUV, EUV, and 
X-ray lines of Fe XVIII and Fe XIX are particularly useful since the strong n=2—*>2 lines 



are essentially entirely produced by direct collisional excitation, and thus should be easier 
to interpret than FUV or X-ray lines, which can include contributions from other processes, 
such as proton impact excitation and dielectronic recombination (DR). 


2. Spectral Models 

We use the Astrophysical Plasma Emission Code version 1.3 (APEC, Smith et ah 2001) 
to predict the Capella spectrum.^ The APEC models for Ee XVIII and Fe XIX contain 501 
and 994 hne-structure levels, respectively, up to principal quantum number n = 5. They 
include the effective collision strengths and atomic transition probabilities calculated using 
the Hebrew University Lawrence Livermore Atomic Code (HULLAC, Liedahl et al. 1995). 
For Fe XVIII, the collision strengths for the ^^* 1/2 ^ ^^ 3/2 transition include resonance 

excitation from R-matrix calculations (Herrington et al. 1998). Proton impact excitation 
rates within the ground state are included for Fe XVIII (Foster et al. 1994) and Fe XIX 
(R. Reid, private communication, 1999). Laboratory X-ray wavelengths (Brown et al. 2002) 
have been incorporated. APEC currently includes DR rates to excited levels of Fe XVII and 
H- and He-like ions, but not for the other Fe L-shell ions. Similarly, DR satellite lines are 
present in APEC for Fe XVII (Safranova et al 2001), but not for Fe XVIII and Fe XIX. 

^APEC VI.3 models, calculated at the low density limit (W = 1.0 cm“^), and the atomic rate data used 
to produce them are available at http://cxc.harvard.edu/atomdb/. 
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3. Data Analysis 

3.1. Observations and Data Reduction 

Multiple spectra of Capella, acquired between 1999 August and 2002 October include 
pointings with the Chandra High Energy Transmission Grating (HETG) with the AGIS-S 
detector for a total exposure time of 182.2 ks, and with the Low Energy Transmission Grating 
(LETG) and HRG-S detector for a total exposure time of 234.2 ks. The HETG and LETG 
data, obtained from the Chandra archive, were reprocessed using GIAO version 3.0^ with 
only minor deviations from the standard pipeline procedures. Effective areas were generated 
for each dataset using the Chandra calibration database GALDB 2.8 and were exposure-time 
weighted to create average effective areas for the summed spectra. 

Extreme Ultraviolet Explorer {EUVE) spectra obtained in 1999 September, which are 
nearly simultaneous with a Chandra LETG/HRG-S pointing, were processed using standard 
EUVE Guest Observer software (IRAF). Agreement between LETG and EUVE fluxes for 
the lines discussed in this paper is good to within about 5%, and henceforth LETG fluxes 
will be used. The Ear Ultraviolet Spectroseopic Explorer {EUSE) line fluxes are taken from 
spectra of Young et ah (2001). 


3.2. Modeling and Measurements 

We calculate the global continuum spectra produced by bremsstrahlung, radiative re¬ 
combination continuum, and two-photon emission over the observed Chandra spectral range. 
We then £t the temperature of the continuum model to the line-free regions of the HETG 
spectrum, identified both from the APEG line list and by visual inspection, which yields a 


^http://cxc. harvard.edu/ciao/ 



6 


temperature of 6 MK, near the peak of the EMD. Since the LETG spectrum is contaminated 
by high order emission, the same continuum model derived from the HETG data is also ap¬ 
plied to the LETG htting. We adopted the abundances of Brickhouse et ah (2000), who 
found no evidence for deviation from the solar abundances of Anders & Grevesse (1989). 
Individual line fluxes from the Chandra spectra were measured using Sherpa (Freeman et ah 
2001) to ht functions approximating the instrumental line prohles. Plus and minus orders 
were ht separately, with the requirement that the line huxes be the same. A narrow range 
for the FWHM was allowed (for HETG, 0.01 - 0.0135 A, and for LETG, 0.045 - 0.06 A), 
standard binning was maintained, and the Gash statistic was applied (Gash 1979). Table 1 
gives the observed huxes for the Fe XVIII and Fe XIX lines with la errors. 


3.3. Model Assumptions 

A continnous EMD (Brickhouse et ah 2000), which is needed to estimate the contribu¬ 
tion of line blends from ions over the entire temperature range, is normalized to the hnx of 
the Fe XVIII A93.92 resonance line and used to predict the line huxes given in Table 1. We 
note that there is only a few percent diherence between the single temperatnre continuum 
model and the EMD for the lines of interest. Since some Fe XIX line emissivities show mod¬ 
est density-sensitivity between the low density limit and densities expected under coronal 
conditions, we have used the APEG code to compute models for a wide range of densities. 
The most ahected line ratio is that of A101.55 to A108.37. At W = 10^° cm“^ the predicted 
ratio is 0.347 (in photon nnits) compared with 0.261 at the standard APEG low density 
limit, in better agreement with the observed hux ratio of 0.328. 

Lack of signihcant variability fnrther supports the assumption that the plasma condi¬ 
tions are stable, as individual lines of Fe XVII, Fe XVIII, and Fe XIX show modest hux 
changes (< 10% deviation from the average value) between Chandra pointings and the 



lightcurves show low levels of variability (<8%) during a single pointing. There is also no 
evidence to challenge the standard assumptions of negligible optical depth (Canizares et ah 
2000; Brown et ah 1998). 


4. Results and Discussion 

Figure 1 compares the observed Fe XVIII and Fe XIX line fluxes to those predicted 
by the spectral codes, APEC, CHIANTI V4.2 (Dere et al. 1997; Young et ah 2003), and 
SPEX VI. 1 (Kaastra et al. 1996), which incorporates the MEKAL model (Mewe et al 
1995). Emissivities provided by M. F. Gu (private communication, 2004) using the Flexible 
Atomic Code (FAC, Gu 2003) are also compared. The fluxes are scaled by the fluxes of their 
respective strong EUV resonance lines, for which direct excitation dominates. All models in 
the hgure are calculated at a single temperature Te = 6 MK, and the same density N^. = 10^° 
cm“^, except for SPEX, which is at the zero density limit. 

Most striking is the discrepancy between the EUV and X-ray lines: the observed X-ray 
fluxes are stronger than predicted fluxes in all models. Even the X-ray 3d-2p resonance 
lines, Fe XVIII A14.208 and Fe XIX A13.518, are under-predicted relative to their EUV 
counterparts by more than 30% and a factor of 2, respectively. Since these factors are larger 
than expected from calibration errors or line blending, it is possible that the accuracy of 
the direct excitation rate coefficients might explain the predicted weakness of AI4.208 (see 
Brown et al. 2005); however, it is difficult to reconcile that with the larger discrepancy for 
AI3.5I8. 

The Fe XVIII and Fe XIX FUV forbidden line fluxes are in good agreement with the 
EUV line fluxes of A93.92 and AI08.37 for the APEC models, and somewhat better than for 
the FAC rates. FAC does not calculate proton impact excitation rates, which are included 



in both APEC and CHIANTI. In APEC models, proton impact excitation increases the 
forbidden line emissivities by 15% and 8% for A974.86 and A1118.07, respectively. The 
predicted FUV line fluxes also begin to increase with density above Nf. 10^2 cm-T APEC 
models give the best agreement at A^^ = 2 x 10^^ cm“^, but are also consistent within 
observational errors with the lower coronal density range. 

Figure 1 also shows some large discrepancies among the strongest X-ray lines, reflecting 
the 3s-2p/3d-2p pattern. For these transitions the largest difference among the predictions 
results from the number of processes calculated with each model. For example, even though 
APEC and CHIANTI have similar collision strengths for the Fe XVIII A15.625 line, additional 
line flux in APEC is produced by direct excitation to n = 4 and 5 levels, followed by radiative 
cascades, while CHIANTI currently includes only levels only up to n = 3. On the other hand, 
for Fe XVIII A16.07 APEC and CHIANTI both show differences of more than a factor of 2 
from FAC because neither includes the effects of DR on the upper level population, which 
are included in FAC. 

Comparisons of APEC and FAC predictions to the observed fluxes of the X-ray lines 
listed in Table 1 are shown in Figure 2. We conhrm a general 3s-2p/3d-2p discrepancy 
pattern for APEC models that is largely removed with the FAC calculations. The 3s-2p/3d- 
2p ratios of the summed line fluxes from APEC are smaller than the observed ratios by 
~20%, whereas FAC agreement is within 10%. The inclusion of DR in the FAC models 
produces the additional 3s-2p line emissivity. 

Another signihcant disagreement between the models and observations occurs for ra¬ 
diative transitions that terminate on excited levels, namely Fe XVIII A15.870, A16.159, and 
A17.623 and Fe XIX A15.198 and A16.110. Although the APEC line list, which is reasonably 
complete in this spectral region, does not include DR satellite lines from either Fe XVIII or 
Fe XIX, blending with satellite lines or lines from other ions cannot explain the extent of the 
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under-prediction. It is possible that the large theoretical wavelength inaccuracies for these 
lines, up to a few percent, have led to misidentihcations in the laboratory measurements. 
Blending of nearby lines from the same ion could produce such a pattern of under and over¬ 
prediction. For Fe XVIII A15.870, this latter explanation is consistent with new wavelength 
calculations (Kotochigova et ah 2005; M. F. Gu 2005). 


5. Conclusions 

A surprising result of this benchmark spectral modeling study is the large discrepancy 
between modern theory and the Capella observations for the X-ray and EUV resonance 
lines of Fe XVIII (30%) and Fe XIX (factor of 2). New FAC calculations including dielec- 
tronic recombination bring most X-ray lines into good agreement with observations; however 
puzzling discrepancies as large as a factor of 2 still remain for some relatively strong lines. 
Additional laboratory and theoretical work is needed to eliminate the largest remaining 
problems. Meanwhile, errors can largely be minimized by judicious choice of line diagnostics 
and consideration of appropriate atomic processes. 

This work is supported in part by the Chandra X-ray Observatory Center (NAS8-39073). 
We thank the CXC staff, particularly Harvey Tananbaum, for supporting efforts to obtain 
these data. We also acknowledge the efforts of the developers of the other public spectral 
modeling codes SPEX and CHIANTI as well as M. F. Gu for the atomic structure code FAC. 
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Fig. 1.— Observed to predicted flux ratios of strong lines in the X-ray, EUV, and FUV 
spectral regions. Shown for comparison are the ratios obtained using the APEC, CHIANTI, 
and SPEX spectral codes and the FAC rates. The density is N^. = 10.0 cm“^, except for 
SPEX. Top: Comparison for Fe XVlll lines, normalized to A93.92. The X-ray lines plotted 
here are A14.208, A15.625, and A16.071. Bottom: Comparison for Fe XIX lines, normalized 
to A108.37. The X-ray lines plotted are A13.518, A14.664, and A15.079. 
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Fig. 2.— Observed to predicted flux ratios of X-ray lines using FAC and APEC. Lines from 
Table 1 excluding heavily blended Fe XVIII A16.004 are shown. Note the 3d-2p lines are 
between 14 and 15 A for Fe XVIII and shortward of 14 A for Fe XIX. Ratios are calculated 
at Ng = 10.0 cm“^. Dash-dotted lines represent agreement within a factor of 2. Top: 
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Table 1. Fe XVIII and Fe XIX Line Measurements 


Inst 

Ion 

^ref 

(A) 

^obs 

(A) 

Transition 

Ju - Jl 

Model Flux^ 

( ph cm“^ ks“^) 

Observed Flux 

( ph cm“^ ks“^) 

fuse'' 

Fe XVIII 

974.86 

974.85 

2p^ 

1 3 

2 2 

5.063 

5.50 ± 0.03 

LETG 

Fe XVIII 

103.93 

103.98 

2s2pe 2 Si/2 - 2p^ 2Pi/2 

1 1 

2 2 

1.625 

1.69 ± 0.05 

LETG 

Fe XVIII 

93.923 

94.02 

2s2p0 2Si/ 2 - Zp'' ^P3/2 

1 3 

2 2 

4.441 

4.44 ± 0.03 

MEG 

Fe XVIII 

17.623 

17.620 

2p^3p 2p3/2 - 2s2p<^ 2Si/2 

3 1 

2 2 

0.300 

0.30 ± 0.01 

MEG 

Fe XVIII 

16.159 

16.163 

2s2p^3s ^P 3/2 - '2s2p^ ^Si/2 

3 1 

2 2 

0.164 

0.13 ± 0.00 

MEG"" 

Fe XVIII 

16.071 

16.073 

2p«(3p)3s «P5/2 - Zp" =P3/2 

5 3 

2 2 

0.418 

0.82 ± 0.01 

HEG"" 

Fe XVIII 

16.071 

16.076 

2p«(3p)3e «P5/2 - Zp" =P3/2 

5 3 

2 2 

0.418 

1.00 ± 0.06 

HEG 

Fe XVIII"' 

16.004 

16.008 

2p*i^P)3s ZP 3/2 . 2p^ ZP 3/2 

3 3 

2 2 

0.768 

0.81 ± 0.04 

HEG 

Fe XVIII 

15.870 

15.873 

2p*Cd)3s ^0^/2 - 2p^ Zpi/2 

3 1 

2 2 

0.095 

0.34 ± 0.02 

HEG 

Fe XVIII 

15.824 

15.831 

2p^(3p)3e «P3/2 - Zp5 ZP 3/2 

3 3 

2 2 

0.179 

0.29 ± 0.02 

HEG 

Fe XVIII 

15.625 

15.628 

2p*Cd)3s ^D^^2 - Zp'' ^^3/2 

5 3 

2 2 

0.290 

0.43 ± 0.02 

HEG 

Fe XVIII 

14.571 

14.559 

2p*(3p)3d '‘P 3/2 - Zp" ZP 3/2 

3 3 

^ 

0.110 

0.21 ± 0.07 

HEG 

Fe XVIII 

14.534 

14.539 

2p*(3p)3d Zpg^j . 2p^ ZP 3/2 

5 3 

2 2 

0.210 

0.39 ± 0.09 

HEG 

Fe XVIII 

14.373 

14.376 

2p'‘(3p)3d ^D^^2 - 2p^ ^P3/2 

5 3 

2 2 

0.278 

0.55 ± 0.02 

HEG 

Fe XVIII® 

14.256 

14.261 

2p‘^CD)3d ZSi/2 - Zp" =P3/2 

1 3 

2 2 

0.087 

0.42 ± 0.03 





‘^Pl/2^P3/2^^b/2 - ^-^3/2 

5 3 

2 2 

0.141 


HEG 

Fe XVIII 

14.208 

14.208 

‘^Pl/2^P3/2^^b/2 - ^-^3/2 

3 3 

2 2 

0.381 

1.40 ± 0.05 





2p^(lD)3ti ''D 5/2 - 2p^ ZP 3/2 

5 3 

2 2 

0.695 


HEG 

Fe XVIII 

11.527 

11.528 

2Pi/2^P3/2'^'^5/2 - ^-^3/2 

5 3 

2 2 

0.032 

0.17 ± 0.01 





2p*fp)4d ZP 5/2 - Zp'' =P3/2 

5 3 

2 2 

0.061 


HEG 

Fe XVIII 

11.423 

11.424 

2p*(3p)4d Zpg^j . 2p5 ZP 3/2 

5 3 

2 2 

0.080 

0.13 ± 0.01 


Fe XXII 

11.427 


2s2p-^/2Sp^/2 - 2p ^Fi/2 

3 1 

2 2 

0.007 


HEG 

Fe XVIII 

11.326 

11.327 

2p'‘(lD)4d ZSi/2 - Zp'' =P3/2 

1 3 

2 2 

0.019 

0.13 ± 0.02 





2p'‘(lD)4d ZP 3/2 . 2 p'' ZP 3/2 

3 3 

2 2 

0.031 






2p'‘(lD)4d ZP 5/2 . 2 p'' ZP 3/2 

5 3 

2 2 

0.038 


FUSE 

Fe XIX'>"^ 

1118.07 


2p* 3 Pi - 2p^ 3P2 

1-2 

1.833 

1.74 ± 0.22 

LETG 

Fe XIX 

120.00 

120.04 

2s2p3 3p2 . 2p^ 3p^ 

2-1 

0.836 

0.97 ± 0.03 

LETG 

Fe XIX 

111.70 

111.74 

2s2p3 3p3 . 2p^ 3p^ 

1-1 

0.326 

0.46 ± 0.02 

LETG 

Fe XIX 

109.97 

109.99 

2s2p3 3 Pi - 2p^ 3p|^ 

1-0 

0.413 

0.46 ± 0.02 

LETG 

Fe XIX 

108.37 

108.39 

2s2p3 3p2 . 2p4 3p2 

2-2 

3.091 

3.13 ± 0.05 

LETG 

Fe XIX 

101.55 

101.59 

2s2p3 3p3 . 2p^ 3p2 

1-2 

0.838 

1.02 ± 0.03 

LETG 

Fe XIX 

91.02 

91.054 

2s2p3 'Pi - 2p^ ^D2 

1-2 

0.241 

0.45 ± 0.02 

HEG 

Fe XIX 

16.110 

16.111 

2pi/2Zp3y23Pl/2 - Zs2p3 ^P 2 

2-2 

0.120 

0.14 ± 0.03 

HEG 

Fe XIX 

15.198 

15.204 

2p\i22p\i2fls - 2s2p3 3p2 

2-2 

0.080 

0.39 ± 0.02 

HEG 

Fe XIX 

15.079 

15.083 

2p3(^S)3s 3S2 . 2p^ 3p2 

2-2 

0.094 

0.33 ± 0.02 

HEG 

Fe XIX 

14.664 

14.671 

2p3(3_D)3a 3p^ . 2p4 3p2 

3-2 

0.079 

0.21 ± 0.01 

HEG 

Fe XIX 

13.795 

13.795 

2Pl/22P3/2^'^5/2 ■ 2^“^ ^ P 2 

3-2 

0.105 

0.24 ± 0.02 





2p3(3_D)3d 3p^ _ 2p4 3p2 

3-2 

0.012 


HEG 

Fe XIX 

13.518 

13.523 

2p3(Z_D)3d 3p^ . 2p4 3p2 

3-2 

0.262 

0.52 ± 0.03 

HEG 

Fe XIX 

13.497 

13.507 

Zpi/2Zp3y23‘*3/2 " Zp^ 3p2 

2-2 

0.118 

0.32 ± 0.02 


Fe XXI 

13.507 

13.507 

lsZ2s2pZ^23s - Is32s2p3 3p)3 

2-2 

0.025 


HEG 

Fe XIX 

13.462 

13.470 

2p3(Z_D)3d 3g^ . 2p4 Sp^ 

1-2 

0.072 

0.25 ± 0.02 

HEG 

Ne IX 

13.447 

13.446 

IsZ 'So - ls2p'Pi 

1-2 

0.397 

0.40 ± 0.02 
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Line blends are listed separately if they contribute > 10% to the Fe line of interest (in the model). Fluxes (including blends) normalized 
to the A93.92 line predicted by the emission measure distribution model using APEC at a density of 1 cm~^ are listed. The observed 
fluxes have been corrected for interstellar absorption using N= 1.7 X 10^® cm~^ (Piskunov 1997), neutral helium, and H/He abundance 
ratio set at 11.6 (Kimble et al. 1993). The largest correction at A120.0 amounts to only 9%. 

^ See Young et al. (2001). 

^ MEG and HEG measurements of this line are given to show the cross-calibration. HEG is preferred for this analysis because of its 
better spectral resolution. 

^ Contribution of OVIII to this line is more than 50%. 

® LETG flux was measured to cross check the calibration of LETG vs HETG. The LETG line is somewhat blended, but the flux is 
within 30% of the HETG flux. 

^ This FUSE measurement is uncertain as this line is blended. Solar network spectra were used to estimate the contribution of C I to 
the blend. 



